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Abstract 

We obtain variational formulas for holomorphic objects on Riemann sur¬ 
faces with respect to arbitrary local coordinates on the moduli space of 
complex structures. These formulas are written in terms of a canonical 
object on the moduli space which corresponds to the pairing between the 
space of quadratic differentials and the tangent space to the moduli space. 
This canonical object satishes certain commutation relations which appear 
to be the same as the ones that emerged in the integrability theory of 
Whitham type hierarchies. Driven by this observation, we develop the the¬ 
ory of Whitham type hierarchies integrable by hydrodynamic reductions 
as a theory of certain differential-geometric objects. As an application we 
prove that the universal Whitham hierarchy is integrable by hydrodynamic 
reductions. 
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1 Introduction 


Various calculations with holomorphic objects on a Riemann surface can be done efficiently 
by using the Fay identity mi. It seems that the Fay identity contains all information about 
identities between the Riemann theta function, normalized holomorphic differentials, the prime 
form and their derivatives with respect to coordinates on a Riemann surface and on its Jacobian. 
On the other hand, these holomorphic objects also depend on moduli of complex structures 
and one needs to be able to compute variations with respect to these moduli. Such formulas 
were obtained by Rauch [3]. He represented a Riemann surface as a ramified covering of CP^ 
and computed variations of holomorphic objects in terms of branch points of this covering. 
Rauch formulas have proven its usefulness and efficiency in various contexts [H m E]. It is 
desirable however to have universal variational formulas which are independent of a partic¬ 
ular representation of a Riemann surface and work for arbitrary coordinates on the moduli 
space. This problem can be approached as follows. It is known that the space of quadratic 
holomorphic differentials on a Riemann surface £ is dual to the tangent space of the moduli 
space Mg of complex structures at the point corresponding to £ (see [7] for a general theory of 
deformations of complex structures and [8] for the Serre duality theorem). Let ui, ...,V 3 g -3 be 
local coordinates on Mg, let be the corresponding basis in the tangent space and 

gi{p)dp^,g 3 gs{p)dp‘^ be the dual basis in the space of quadratic differentials. The object 


33-3 

G{p)dp^ = ^ gi{p)dp^ 

i=l 


A 

dvi 


does not depend on any choice of coordinates. Moreover, let Mg^n be the moduli space of 
Riemann surfaces with n punctures Ui, ...,Un- Here we can vary both the complex structure of 
£ and points ui, ...,«„ in £■ A basis in the tangent space looks like ..., ..., ^ 

and the corresponding object is 


G{p)dp^ = y^^F{p, Ui)dp^ 

i=l 


A 

dui 


^9 3 ^ 

i=i 


where F{p,u) has a pole of order one ai p = u and is holomorphic outside the diagonal. The 
residue of F{p,u) at p = u is a constant and without loss of generality we assume F{p,u) = 
+ 0(1). Roughly speaking, F{p,u) should be a quadratic differential with respect to p 
and a vector field with respect to u. However, the transformation law for this object is more 
complicated. Indeed, if we change coordinates by p = p{p,vi, ...,V 3 g- 3 ), Ui = fi{ui,vi, ...,V 3 g- 3 ) 
(and do not change vi, ...,V 3 g- 3 ), then G{p)dp^ should transform as a vector held with respect 
to Ml, ...,Mji,Mi, ...,M 3 g _3 (and as a quadratic differential with respect to p). This leads to the 
following transformation law for F{p,u): 

A(p,m) = ^^-^(f{p{p),p{u)) -G(p(p))(p(m))) 
p [u) V / 
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with the same coefficient at F as if it were a quadratic differential with respect to the first argu¬ 
ment and a vector field with respect to the second one but with an additional term depending 
on G. 


The first main result of the paper is the set of explicit formulas for the action of G{p) 
on various holomorphic objects on the Riemann surface £ such as prime form, holomorphic 
1-forms, period matrix, see H2.1d\) . ^2.14^ , H2.15\) . 

We have also obtained commutation relations for G{p): 

[G(Pi),G(P2)1 = 


= ■F(P2,Pi)G'(pi) - -F(Pi,P2)G'(p2) + 2F(p2,Pi)„G(pi) - 2F(pi,p2)„G(p2). (1.1) 

The same relations hold for G{p). One can check that these relations are invariant with respect 
to an arbitrary change of coordinates by virtue of transformation laws of G and F. 

Our main motivation for these studies came from attempts to understand better integrable 
structures of the so-called Whitham type hierarchies [9i[ini[n]. Recall that a Whitham type 
hierarchy is defined as compatibility conditions of the following system of PDEs: 

dijj 

— = hi{z,ui,...,Un), i = l,...,N. (1.2) 

Here ...,Un are functions of times ti,...,tN and 2 ; is a parameter. The system (11.21) is 

understood as a parametric way of defining iV — 1 relations between partial derivatives i = 
1, ...,N obtained by excluding z from these equations. Functions hi are called potentials of this 
Whitham type hierarchy. 

An important class of such hierarchies associated with the moduli space of Riemann surfaces 
of genus g with n punctures (the so-called universal Whitham hierarchy) was constructed and 
studied in mm The universal Whitham hierarchy is important in the theory of Frobenius 
manifolds [13], matrix models and other areas of mathematics. Note that the set of times in the 
universal Whitham hierarchy coincides with a set of meromorphic differentials on a Riemann 
surface (holomorphic outside punctures), and that the potentials hfiz) are integrals of these 
differentials. 

A natural question is in which sense a Whitham type hierarchy is integrable. In this paper 
we concentrate on an approach to integrability theory of such systems based on the so-called hy¬ 
drodynamic reductions [TTl [151CB] • In this approach a quasi-linear system is called integrable if 
it possesses a large family of hydrodynamic reductions. This family of hydrodynamic reductions 
must be parametrized by solutions of another system of PDFs called Gibbons-Tsarev system 
[niilHlilg]. Therefore, Gibbons-Tsarev systems play a crucial role in this integrability theory. 
See [T9| and references therein for the definition and examples of Gibbons-Tsarev systems. 

The second main result of the paper is a description of integrable structures that appeared 
in the hydrodynamic reduction method as a certain differential-geometric structure. We call it 
a GT structure. 
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By definition, a GT strnctnre is defined locally by a family of vector fields g{p) and a 
function f{pi,P 2 ) satisfying relations similar to fll.ip . see Section 3 for precise definitions. It 
becomes transparent from these definitions that a natural GT structure exists on the moduli 
space Mg^ri and is represented by the objects G(j)) and F{p,u) described above. 

Given a GT structure one could ask how to hnd all corresponding integrable Whitham type 
hierarchies. It turns out that in order to classify all possible integrable hierarchies with given 
GT structure one needs to hnd all functions X{pi,P 2 ) satisfying the functional equation 

9{Pi){HP2,P3)) = HPuP3)HP2,Pi)pi -HP2,P3)f{PuP2)p2- 
-f{PuP2)KP2,P3)p2 - f{Pl,P3)HP2,P3)p3- 

Moreover, to hnd all potentials h{p) of a given hierarchy one needs to solve another functional 
equation 

9{Pi){h{P2)) = HPi,P2)h'{pi) - f{puP2)h'{p2). 

It is natural to ask if the universal Whitham hierarchy is integrable by hydrodynamic re¬ 
ductions. One could expect that the corresponding GT structure is given by G{p) and F{pi,p 2 ) 
and needs to hnd a function X{pi,P 2 ) which gives the universal Whitham hierarchy. 

The third main result of this paper is a proof that the universal Whitham hierarchy is in¬ 
deed integrable in all genera by hydrodynamic reductions. We give the corresponding function 
A(pi,P 2 ) and the precise form of potentials. 

Let us describe the content of the paper. In Section 2 we recall main dehnitions and notations 
of holomorphic objects on a Riemann surface, construct our main object G{p) and compute its 
action on holomorphic objects. We compute commutation relations for G{p) as well. We give 
some examples and explain how the Rauch formulas are connected with ours. In Section 3 we 
introduce GT structures and develop a theory of these structures. In particular, we explain how 
to construct new GT structures from a given one and how to construct potentials if we are given 
a function \{pi,P 2 ) dehning our hierarchy. We also explain a relation between GT structures 
and Lie algebroids of a certain type. In Section 4 we recall the dehnition and basic properties 
of Whitham type hierarchies. In Section 5 we recall dehnition of Gibbon-Tsarev systems and 
prove that there exists a one-to-one correspondence between Gibbons-Tsarev systems and GT 
structures. In Section 6 we discuss the dehnition of integrability of Whitham type hierarchies in 
our framework. We explain that the dehnition based on hydrodynamic reductions and Gibbons- 
Tsarev systems is equivalent to ours. We refer to [19] for a full discussion of integrability of 
Whitham type hierarchies based on hydrodynamic reductions and Gibbons-Tsarev systems. 
It was not possible to make this paper self-contained and repeat this discussion here without 
essential increasing of the length of the present paper. In Section 7 we recall the dehnition of 
the universal Whitham hierarchy and prove that this hierarchy is integrable by hydrodynamic 
reductions. 
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2 Holomorphic objects on Riemann surfaces and defor¬ 
mations of complex structures 


Let = D/r be a compact Riemann surface of genus 5 ^ > 1, D C C its universal covering and 
r = Tii{£). Denote by Oq,, 6 q,, a = 1 ,(7 a canonical basis in the homology group Hi{£, Z). Let 
us choose a coordinate in D and use the same symbols for holomorphic objects on S and their 
lifting on D. We will also use the same symbol for a point in S, its lifting in D and its coordinate. 
Let uja{z)dz be the basis of holomorphic 1-forms on £ normalized by oj^dz = daji- Choose 
a basepoint zq and dehne the Abel map qa{z) = uja(z)dz. Note that cZa = Qa- Denote the 
prime forrr0 by B(x, q)(dx)~^^^(di/)~^^^. Let updz be the matrix of 6 -periods. Details 

on holomorphic objects on Riemann surfaces are given in [DEI [5]. Recall that 

E{v,u) =—E{u,v), E{u,v) = u — V — ^S{u){u — v)^ + 0{{u — vY), (2.3) 

where S{p) is the Bergman projective connection on £. Note that E{u,v) is multivalued. If u 
or V is moved by a^, it remains invariant. If u moves by to h or n moves by to h, then 


E{u, v) = E{u, v) exp ^ - JiiBaa + 27ri(ga(n) - qaiu))^, 

E{u, v) = E{u, v) exp ^ — mBaa — ‘^^miqai.v) — qa{u))^. 

Let W{u,v) = (\n{E{u,v))uv be the Bergman kernel. Recall that 

/ W{u,v)du = 0, / W{u,v)du = 2'Kiuja{v), / / W{u,v)dudv = 27riBai3. 

J dj J b(y J bct J i) ft 


( 2 . 4 ) 


( 2 . 5 ) 


Recall a description of the tangent space to the moduli space Mg of Riemann surfaces at the 
point corresponding to £ [201 ED- Let p ^ £ he the center of a small disc D <Z £. Let L be the 
Lie algebra of holomorphic vector helds on D \ {p} and Lp, Lout be subalgebras of L consisting 
of vector helds holomorphic at p and holomorphic on \ {p} correspondingly. It is known that 
the tangent space to the moduli space Mg is isomorphic to the quotient Lj (Lp©Lo„t). Let Mg^i 
be the moduli space of Riemann surfaces with a puncture at u E £. The tangent space to Mg 1 
is isomorphic to the quotient L/{Lp © Lout,u) where Lout,u C Lout consists of vector helds with 
zero at u. Let us construct vector spaces dual to these tangent spaces using the Serre duality 
theorem [S]. There exists a non degenerate pairing between the space L and the space Q of 
quadratic diherentials holomorphic on D \ {p}. This pairing is given by {v, q) = ReSp(ng). The 
space dual to the tangent space of Mg is equal to (Lp © Lout)'^ C Q and consists of quadratic 
diherentials holomorphic on £. Similarly, the space dual to the tangent space of Mp 1 is equal 

^In this paper we represent differential-geometric objects as functions with prescribed transformation laws 
with respect to an arbitrary change of coordinates. For example if a; = f{x), y = then the prime form 

transforms as E{x,y) = y'{x)~^^'^y,'{y)~^^'^E{y{x), fi{y)). 
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to {Lp © Lout,u)^ C Q and consists of quadratic differentials holoniorphic on \ {«} with pole 
of order less or equal to one at u. More generally, the space dual to the tangent space of 
Mg^n of the moduli space of Riemann surfaces with punctures at ui, ...,Un consists of quadratic 
differentials holoniorphic on \ ...,Un} with poles of order less or equal to one at mi, Un- 

Let ni,...,n 3 g _3 be local coordinates on moduli space Mg. Let ^ be the corre¬ 

sponding basis in the tangent space and gi{p)dp^,g^g_^{p)dp^ be the dual basis in the space 
of quadratic differentials. The object^ 


3g-3 

G{p)dp^ = ^ gi{p)dp^ 

i=l 


A 

dvi 


does not depend on the choice of coordinates. A similar construction for Mg^n gives the object 


G{p)dp^ = '^F{p,Ui)dp^ 

i=l 


A 

dui 


3g-3 

+ '^gj{p)dp^ 

i=i 



where ui, ...,Un are coordinates of n points in S and 

F{pi,P2) = —^— + 0 ( 1 ). 

Pi -P 2 


( 2 . 6 ) 


Proposition 2.1. Under an arbitrary change of coordinates of the form 

p =/i(p,ni, ...,n3g_3), Ui = iJ,{ui,Vi,...,V3g_3) (2.7) 

the objects G{p), F{pi,p 2 ) obey the following transformation rules 


G{p)=p'{prG{p{p)), ( 2 . 8 ) 

Fipi,P2) = ^^^-(F(p(pi),/i(p2)) - 0(p(pi))(/i(p2))) (2.9) 

Proof. The relation fl2.8l) means that G{p) is a quadratic differential in p (with values 
in vector helds in vi, ...jV^g^^). In order to obtain fl2.9l) we perform an arbitrary change of 
coordinates of the form p = /i(p, ui, ..., ^ 33 - 3 ), Ui = p,{ui, vi, ..., v^g-^), Vj = Vj and require that 

the object G{p) transforms as a vector field in ui, ...,Un,vi, ...,V3g-3. The relation (12.91) is a 

consequence of this requirement. □ 

Proposition 2.2. The following identities hold 


[G(pi).G(k)] = -F(p 2 ,Pi)G'(p,) - -F(p,.P 2 )G'(p 2 )+ (2.10) 

^Note that the functions G, gt, F etc. depend also on vi, ...,1133-3. We will often omit these arguments in 
order to simplify formulas. 
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+2F{p2,pi)p,G{pi) - 2F{pi,p2)p^G{p2), 


[Gipi),G{p2)] = F{p2,p,)G\p,) - Fip,,p2)G\p2)+ ( 2 . 11 ) 

+2F{p2,pi)p,G{pi) - 2F{pi,p2)p^G{p2), 


G{P2){F{PI,P3)) - G{pi){F{p2,P3)) = F{pi,p2)F{p2,P3)p2 - F{p2,Pi)F{pi,p3)p^+ (2.12) 

+ F{puP3)F{p2,P3)p3 - F{p2,P3)F{pi,p3)p^ + 2F{p2,P3)F{pi,p2)p, - 2F{pi, p3)F{p2, Pi)p„ 


G{pi){E{p2,P3)) 

E{P2,P3) 


^F{pi,P2)p^ + ^F{pi,P3)p^- 


-F{pi,P2) 


E{P2,P3)p2 

E{P2:P3) 


F{pi,P3) 


E{P2,P3)p3 

E{P2,P3) 


1 / E(Pi,P2)pi 

2 V E{pi,p 2 ) 


E{j)l,P3)p3 y 

E{pi,P3) ) 


(2.13) 


G(pl) 



F{px,P‘^UJi{p2) - F{px,P3)^i{p3)- 


F{p\,P‘^P3 

F{pi,P2) 


^i{pl) + 


E{J)1,P3)p3 

E{puP3) 


^i(pl), 


(2.14) 


G{p){Bjk) = 2mujj{p)uJk{j)). (2.15) 

Proof. Notice that (12.111) is a formal consequence of (I2.10j) and (I2.12p (see Proposition 
3.1). 

Consider the difference of the l.h.s. and the r.h.s. of each of fl2.10p . fl2.12p . fl2.13p . fl2.14p . 
Expanding these expressions on each diagonal Pi = Pj, i ^ j and using fl2.3l) and (12.61) one 
can check that each of these expressions is holomorphic on all diagonals. Making an arbitrary 
change of coordinates of the form pi = /i(pi, Ui,..., fsg-s), z = 1, 2, 3 one can check that all these 
differences are transformed as tensor fields in pi,P 2 ,P 3 - In particular, the difference between the 
l.h.s. and the r.h.s. of fl2.12p is a holomorphic quadratic differential in pi,p 2 and holomorphic 
vector field in P 3 . This proves fl2.12p because any holomorphic vector field vanishes. Similarly, 
the differences between the l.h.s. and the r.h.s. of (I2.13p . (I2.14p are holomorphic quadratic 
differentials in pi and holomorphic functions in P 2 ,P 3 - Moreover, these functions vanish on the 
diagonal p 2 = P 3 - This would prove (12.131) . (I2.14p (any holomorphic function is a constant) 
provided that we prove that the differences between the l.h.s. and the r.h.s. are single valued. 

Taking the second derivative of the equation (12.131) we get 


G'(pi)(H'’(P2,Ps)) 


(2.16) 








































P2P3 




E{P2,P3)p 3 E{pi,P2)p,E{pi,P3) 


/Pl 


E{P2,P3) ’ E{p2,P3) E{PI,P2)E{PI,P3) 

where (W{p 2 ,P 3 ) = (ln(£'(p 2 ,P 3 ))p 2 P 3 Bergman kernel. Let us prove this identity. Let 

A{pi,P 2 ,P 3 ) be the difference of the l.h.s. and the r.h.s. of fl2.16p . It is a quadratic differential 
in Pl and 1-form in both P 2 ,P 3 - Using transformation properties (12.dh we see that A{pi,p 2 ,p 3 ) 
is single valued. Therefore, A(pi,p 2 ,P 3 ) = Yl^,i3=i''"ap{Pi)oOa{P2)(^i3{P3) where r^isipi) are some 
holomorphic quadratic differentials. Computing A(pi,p 2 ,P 3 )dp 2 dp 3 we obtain Ta/sipi) = 

0 which proves (12.161) . Computing dp 2 dp 3 of the l.h.s. and the r.h.s. of (I2.16p and using 

fl2.5l) and fl2.4l) we obtain fl2.15l) . The difference between the l.h.s. and the r.h.s. of fl2.14p is 
single valued by virtue of fl2.15l) . This proves fl2.14p . Equation fl2.13p is proven in a similar way. 
Note that the difference between the l.h.s. and the r.h.s. of fl2.13p is single valued by virtue 
of (12.141) . Equation (I2.10p is proven by applying its l.h.s. and the r.h.s. to Bjk. For example, 
on the l.h.s. we have G(pi)(G(p 2 )(-Bjfc)) — G{j) 2 ){G{j)i){Bjk))- Computing by virtue of fl2.15p . 
fl2.14p we prove fl2.10p . □ 

Remark 2.1. Recall that the Riemann theta-function is defined by 
6{zi ,..., Zg) = exp(27rim • z -|- vrimBm*). 


Here we use bold symbols for the corresponding vectors: m = {mi, z = {zi,...,Zg), 

m ■ z = miZi mgZg, and B is the period matrix. We have 


G{p){e{zi,...,Zg))= 

a, 13=1 


de{zi, ...,Zg 


dB. 


a/3 




d'^e{Zl, ...,Zg 


a,13=1 


dZadZ/s 


-Ua{p)ujfl{p) 


where we used heat equation for 9 and fl2.15p . 

Remark 2.2. Expanding (I2.13p on diagonal p 2 = P 3 we obtain 

G{pi){S{p2)) + E{pi,P2)pI + 2S{p2)E{pi,p2)p^ + S{p2)p^E{pi,p2) - QW{pi,P2f = 0. 


Example 2.1. Let g = 2. Represent as a 2-fold covering of CP^. Let x be an affine 
coordinate in CP^ and let branch points of the covering be at x = 0,1, oo, a, b, c. The curve S 
is given by = x{x — l)(x — a){x — h){x — c). One can check that 


G{p) = 


a[a 


2p{p — 1)\ p — a da 


1 ) a ^ b{b 


P 


1 ) d c(c 


1 ) d 


b db 


P 


c dc 


F{pi,P2) = 


{pi - a){pi - b){pi - c)p 2 {p 2 - 1) + qiq 2 


2(pi - P 2 )Pi{Pi - l)(Pi - o){pi - b){pi - c) 
where p,Pi,P 2 are affine coordinates in CP^ and qf = Pi{pi — l){pi — a){pi — b){pi — c). 
Let us compare our variational formulas with Rauch ones. The equation fl2.15p reads 
1 


a(a — 1) d b{b — l) d c(c — 1) d 

2p{p — 1) \ p — a da p — b db p — c dc 


{Bjk)dp^ = 2Tiiujj{p)dp ■ ujk{p)dp. (2.17) 
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Let r be a local coordinate near branch point a, we have p = a + dp = 2rdr. Expanding 
the l.h.s. of fl2.17l) we get (2^^ + 0{T))dT‘^. Therefore = ni\p=a\p=a and we 
arrive at a Rauch formula. 

In general, if S is represented as a branched covering of CP^ ramihed at ak G CP^ with 
ramihcation indexes Vk, k = 1,2,..., then G{p)dp^ = ^ + o((p —and Rauch 

formulas can be derived from ours in a similar way. 

Example 2.2. Let us choose ..., Pjg^ g ^ 3 ^ 3 as local coordinates in Mg. Applying 

(I2.15P to Bj^^ki we get 

3g-3 Q 

G{p) = 27ii ujj^{p)ujkXp)-Q^—- 

Applying again fl2.15p to an arbitrary Bjk we obtain quadratic relations between normalized 
differentials. Namely, if S{Bii, B 12 , ...,Bgg) = 0 is a relation between the entries of the period 
matrix (there are _ 3 ^ _|_ 3 functionally independent ones), then 

j^k = l 

See [ 22 ] for a similar formula for quadratic relations between normalized differentials. 


3 GT structures 

Based on the identities (I2.10p . (12.121) we want to introduce a general differential-geometric 
structure: a family of vector fields g{p) and a function f{pi,P 2 ) satisfying the same relations. 
We will see later that this structure is equivalent to an integrability structure of Whitham type 
hierarchies, the so-called Gibbons-Tsarev system. 

Let g{p) = 5'i(p, Wi) be a family of vector helds parameterized by p and 

f{Pi,P2,Vi,...,Vm) be a function. 

Definition 3.1. A local GT structure is a family g{p) and a function f{pi,P 2 ) satisfying 
the following relations: 


[9iPl),9iP2)] = fiP2,Pi)9'iPi) - /(Pl,P2)^'(P2 ) +2/(p2,Pi)pi^(Pi) -2/(pi,P2)p2^(P2), (3.18) 


9{P2){f{pi,P3)) - 9{Pl){f{P2,P3)) = f{Pl,P2)f{p2,P3)p2 “ f {P2, Pi) f {Pl, P3)pi+ (3.19) 

+ f{PuP3)f{P2,P3)p3 - f{P2,P3)f{PuP3)p3 + 2/(^2, P3)/(Pl, ^ 2)^2 “ 2/(pi, P3)/(P2, Pl)pi, 
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f{Pi,P2) = —-— + 0(1). (3.20) 

P1-P2 

Here and in the sequel we often omit additional arguments Vi, ...Vm, indexes stand for partial 
derivatives and g'{p) = ^ 

Given a GT structure we can construct new GT structures in different ways. 

Proposition 3.1. Let g{p), fiPi,P2) satisfy relations fid.isp . fld.lOp and 

d d 

9{p) = + ■■■ + fip^'^n)^ + g{p)- (3.21) 

Then g{p)^ fiPi,P 2 ) also satisfy relations fl3.18p . (13.1 Op . 

Proof. Equation fl3.18p is verihed by direct computation for n = 1 and through induction 
by n for n > 1. Equation fl3.19p remains the same because f{pi,P 2 ) does not depend on 

Ui,...,Un. □ 

We say that a GT structure given by g{p), f{pi,P2) is obtained from a GT structure g{p), 
f{pi,P 2 ) by adding n points ui, ...,Un- This procedure corresponds to a regular helds extension 
of a Gibbons-Tsarev system HI. 

Proposition 3.2. Let g{p), fipi,P2) satisfy relations (I3.18p . (I3.19p and 


= E 

) * * * 5 
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(3.22) 


{ij,i + + ••• + f{p, Ujfl) d 






du,i, 


+ 9{p)- 




Then fiPi,P 2 ) also satisfy relations (13.181) . (I3.19p . 

Proof. Let us start with the following local GT structure 


d 


9{P)= + ^(p)- 






We make the following change of coordinates 


(3.23) 


'^j,0 ^j,0i 

Vj,2 = Ujfl + 2eMj,i + 


(3.24) 
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nAuj — 1)2 n ■ 

'^j,nj — Ujfi + TljeUj^i H---e Uj ^2 + ••• + e ^Uj^ny 


In new coordinates we have 


9 {p)= y] (fiP^Uyo)- -+-(/(p,Mj,o +-/(p,Mi,o))^-d 

i<i<k ^ e duyi 


1 d 

+ ^(/(p, Wi,o + 2eM«,i + e^Mj-2) - 2f{p,Ujfl + ew^-i) + f{p,Ujfl))- -h ... 


1 d \ 

+ —(/(P, Uj^o + rijeuj^i + ... + e^^Vj^n^) - ... + {-lpf{p, j + ^(p)- 


dUj^2 
d 

j,nj 


Taking the limit e ^ 0 we obtain fl3.22l) . □ 

We say that the GT structure (13.221) is obtained from the GT structure fl3.23p by colliding 
points Vjfi, Vj^i ,..., Vj^nj for each j. 

Remark 3.1. Equation (13.191) is equivalent to Jacobi identity for (I3.18p provided that vector 
£eldsfi((pi), g{p2), gipa), g'ipi), g'{P2), 9 '{Pz) are linearly independent for generic pi, pa, Pa- 

Remark 3.2. A local GT structure can be regarded as a certain Lie algebroid. Let 


g{p) = 62 + (p - z)e2. + (p - zfe^, + 

In other words, let 6^+2 = d g^^\z). Let Ci = ^ and 

1 


/(Pi,P2) = —^ fi,j{z) (pi - zy{p2 - zy 


P1-P2 


i,j=0 


Then we have [ei, e*] = {i — l)ei+i and equation (I3.18P is equivalent to 


i—l j—1 

[G, ej] = (j - OG+j + + r - l)/j_2,rG-r+l - ^(j + r - l)fi-2yej_r+l- 

r=0 r=0 

In particular, if /(pi,P 2 ) = 1 then we get [cj, Cj] = [j — i)eiyj for ei, 62 ,.... Note that (I3.19P 

always holds for /(pi,P 2 ) = p^\p^ ■ Therefore, a local GT structure can be regarded as a certain 
deformation of a Lie algebra with basis Ci, 62 ,... and bracket [e*, Cj] = (j — in the class of 
Lie algebroids. 

Given a local GT structure one wants to classify all Whitham type hierarchies that are 
integrable by hydrodynamic reductions and that correspond to a given Gibbons-Tsarev system. 
It turns out that in order to do this one needs to hnd all functions A(pi,p 2 , fi,..., fm) satisfying 
a certain condition. This can be formalized in the following way: 
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Definition 3.2. An enhanced local GT structure is a family of vector fields g{p), a function 
f{pi,P 2 ) and an additional function \{pi,p 2 ,Vi, satisfying the relations fl3.18p . fl3.19p . 

and 


9{Pi){KP2,P3)) = KPi,P3)KP2,Pi)pi - HP2,P3)f{Pi,P2)p2- (3-25) 

-/(Pi,P2)A(P2,P3)p2 - f{PuP3)HP2,P3)p3, 

A(pi.P2) = — !— + 0(1). 

Pi -P2 

Given an enhanced local GT structure one wants to hnd a vector space of all potentials of 
the corresponding Whitham type hierarchy. In all known examples these spaces are spaces of 
solutions of linear systems of PDEs. However, in the general case we can dehne this vector 
space as a space of solutions of a linear functional equation. 

Definition 3.3. Given an enhanced local GT structure we dehne the corresponding vector 
space of potentials as the space of solutions of the following functional equation for a function 

h{p, Vi, ...,Vm)- 


9{Pi){KP 2)) = KPuP2)h'{pi) - f{pi,P2)h'{p2). (3.26) 

Note that expanding fl3.26p near diagonal p 2 = pi we obtain for h{p,vi, ...,Vm) a system of 
linear PDEs equivalent to fl3.26p . 

The following procedure gives a standard way to obtain solutions of fl3.26p : 

Proposition 3.3. Let 7 be a path in C such that 9{x{t,p2)f{pi,t )) ^ Then 

h{j>) = / \(t,p)dt 

J'f 

is a solution of fl3.26p . 

Proof. Substitute this expression for h{p) into (I3.26p and use (13.251) . Direct computation 
shows that the difference between the r.h.s and the l.h.s. of fl3.26l) is 9{x{t,p2)f{pi,t)) ^^^ ^ 

Let us promote local GT structures to differential-geometric ones. 

Proposition 3.4. Relations fl3.18p . fl3.19p . fl3.20p are invariant with respect to arbitrary 
transformations of the form 


Pi = p{pi, Vi ,..., Vm), 9{p) = p'{pf9{p{p)), 
fiPl,P2) = ^^^^(/(/i(pi),/i(p2)) -^(/X(pi))(/i(p2 ))). 


(3.27) 
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Let TT : M —)■ 5 be a bundle with m dimensional fiber F and one dimensional base B. 

Definition 3.4. A GT structure on vr is a local GT structure on each trivialization for each 
U C B such that for different trivializations these local GT structures are connected by fl3.27l) . 
Here vi, ...,Vm stands for coordinates on F and p is a coordinate on B. 

Proposition 3.5. Relations fl3.25p are invariant with respect to an arbitrary transforma¬ 
tions of the form fl3.27p provided that A is transformed as 


HPup2) =/^'(Pi)A(/i(pi),/i(p2)) 


(3.28) 


Definition 3.5. An enhanced GT structure on vr is an enhanced local GT structure on 
each trivialization for each U <Z B such that for different trivializations these enhanced local 
GT structures are connected by (13.271) . fl3.28p . 

Example 3.1. It is clear from (12.101) . (12.121) that g{p) = G{p), f{pi,P 2 ) = F{pi,p 2 ) is a 
GT structure on the bundle Mg^i —>■ Mg. 

Similar GT structures exist for = 0,1- In the case 5 ^ = 0 we consider the moduli space 
M),n +3 of complex structures on CP^ with punctures in n-|- 3 points. We fix 3 points at 0,1, cxo 
and move other points. The formulas for the corresponding GT structures read 


/(Pl,P2) 


P2{P2 - 1 ) 
{Pl-P2)Pl{Pl - 1)’ 


dip) = 

i=l 


Uiiuj - 1 ) d 
{p - Ui)p{p - 1 ) dui' 


(3.29) 


In the case 5 ^ = 1 we consider the moduli space of complex structures on an elliptic 

curve with punctures in u -|- 1 points. We fix one point at 0 and move other points. We also 
deform the complex structure on our elliptic curve. The space of complex structures is one 
dimensional in this case. We use the modular parameter r with Imr > 0 as a coordinate on 
the moduli space of elliptic curves. The formulas for the corresponding GT structures read 

d ^ d 

fiPi,P2) = p{pi-P2,r) - p{pi), g{p) = + ^(p(p-Mj,r) ( 3 - 30 ) 

j=i J 

where p{p,t) = ^ln( 6 '(p, r)) and 9{p,t) = 

Remark 3.3. In these GT structures we can also collide points and obtain new GT struc¬ 
tures. Moreover, in the case g = 0 (resp. 5^ = 1) we can collide points with 0,1, 00 (resp. with 
0) by doing a substitution similar to (I3.24p . In the case 5 ^ = 0 we can also make an arbitrary 
fractional linear transformation with constant coefficients sending 0 , 1 , 00 to a, b, c and collide 
some of a, 6 , c. 

Remark 3.4. Consider an enhanced local GT structure with g{p) given by (I3.23p . Colliding 
points ...,Vj^nj by substitution (I3.24p and taking the limit e —)■ 0 we can do the same 

substitution and limit in the function A and obtain a new enhanced local GT structure. 
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4 Whitham type hierarchies 


Given a set of independent variables ti, ...,tN called times, a set of dependent variables vi, 
called fields and a set of fnnctions hi{z,vi, ...,Vm), i = called potentials we define a 

Whitham type hierarchy as compatibility conditions of the following system of PDEs: 

dib , , , , 

— = hi{z, ni, * = 1, N. (4.31) 

Here ip,Vi, ...,Vm are fnnctions of times and 2 ; is a parameter. The system fl4.3ip 

is nnderstood as a parametric way of defining iV — 1 relations between partial derivatives 
i = 1,N obtained by eliminating from these eqnations. Let ns assnme that the system 
(14.31 p is compatible. Compatibility conditions can be written as 


E(( 

1=1 


fdhi dhj 

dhj dhi ^ 

\ 1 

(dhj dhk 

dhk dhj'' 


/dhk dhi 

dhi dhk\ dvi\ 

V dz dvi 

dz dvi j 


\ dz dvi 

dz dvi f 

1 

' dti 

V dz dvi 

dz dvi ) dtj) 


(4.32) 

where i, j, k = are pairwise distinct. Let Vi^^k be the linear space of fnnctions in 2 ; 

«nnnnprl W _ 2hi§hi 2hi§hk. _ _ dludh^ 7 = 1 m 

^ oz ovi oz ovi ’ az ovi oz ovi ’ oz ovi oz ovi ’ ’ ’ 

Proposition 4.1. Let be finite dimensional and diml^ j^ = D. Then (14.321) is 
eqnivalent to a hydrodynamic type system of D linearly independent eqnations of the form 


m 


1=1 


(arl{Vi,...,V„i)^ + brl{Vi,...,Vm)^ + Crl{Vi,...,Vm)^] =0, r = l,...,D. (4.33) 

\ ut.i ut^i ut}^/ 


Proof. Let {S'i(z),..., S^iz)} be a basis in Vij^k and 

dh±db± _ db±§h± — _ df^d^ _ q dhj^dlH _ dl^dhj^ _ , q 

dz dvi dz dvi — l^r=l^rlOr, — l^r=l^rlOr, 

Snbstitnting these expressions into (14.321) and eqnating to zero coefficients at Si, ...,Sd we 
obtain (14.331) . □ 

Remark 4.1. In the theory of integrable systems of hydrodynamic type the system (I4.3ip 
is often referred to as a psendo-potential representation of the system (14.331) . 

Remark 4.2. In all known examples of integrable Whitham type hierarchies we have 
n < D < 2n — 1. Therefore, this ineqnality can be regarded as a criterion of integrability. 
However, in this paper we explore another criterion of integrability given by the so-called 
hydrodynamic redaction method. 


5 Gibbons—Tsarev systems 

Gibbons-Tsarev systems are the main ingredient of the approach to integrability of Whitham 
type hierarchies and, more generally, to integrability of qnasi-linear systems of the form (I4.33P 
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based on hydrodynamic reductions. In this approach hydrodynamic reductions of a given 
hierarchy are parameterized by solutions of a Gibbons-Tsarev system. In this Section we 
explain a connection between Gibbons-Tsarev systems and GT structures. 

Let pi, ..,Pm, Vi, ■■■,Vm be functions of auxiliary variables ri, and di = 

Definition 5.1. A Gibbons-Tsarev system is a compatible system of partial differential 
equations of the form. 

diPj = f{pi,pj,vi,...,Vm)diVi, = 

diVj = gj{pi,vi,...,Vm)diVi, j = 2,...,m, i = (5.34) 

didjVi = q{pi,Pj,vi,Vm)diVidjVi, i ^ j, i,j = 1,... ,M. 


Remark 5.1. It follows from the compatibility assumption that the space of solutions of 
a Gibbons-Tsarev system is locally parameterized by 2M functions in one variable. Note that 
/, Qi, q do not depend on M and therefore M can be arbitrary large for a given Gibbons-Tsarev 
system. 

We say that a Gibbons-Tsarev system is non-degenerate if f{pi,P 2 ,Vi, has a pole of 

order one on the diagonal p2 = Pi- Assume in the sequel that all Gibbons-Tsarev systems are 
non-degenerate. 

Proposition 5.1. There exists a one-to-one correspondence between non-degenerate Gibbons- 
Tsarev systems and local GT structures. 

Proof. Redehning f,gi from (15.341) we write a Gibbons-Tsarev system in the form 


o / ipi, Pj,Vl, ■ ■ ■ , Um) o ■ _L ' ■ ' 1 yi/f 

OiPi =-- ^OiVx, It 3, b J = 1, • • •, 

gx[j)i,Vx, ...,Vm) 

diVx _ diVj 




gi{pi,Vi, gj{Pi,Vu 

didjVi = q{pi,Pj, vi,..., Vm)diVidjVi, ij^j, i,j = 1,... ,M 


(5.35) 


where f{pi,P2) = p^-p^ + ^(l)- Indeed, is the residue of f{pi,Pj) from fl5.34p at pj = pi 
Write 

^ d 

dip) = ^9iiP,Vl,...,Vm) - 

i=l 


dvi' 


Compatibility of the system fl5.35l) implies c?icl20(P3, ..., w^) = <92^i0(P3, ..., n^) for an 

arbitrary function <f>. This can be written as 

/(P..P2) A + /(P..P3) A + ((/(P3.P3) A + 9(p3))^(P3) . ^) ■ = 


d o \ / o 

fiP2,Pi)^ + fiP2,Pz)^ + 9ipi)j (ifipi,pz)^ + gipMiP'i) ■ 


d 


d 


diui \ d2Ui 


dpi 


dpz 


dpz 


QiiPi)^ 9iiP2)' 
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Expanding this eqnation and equating coefficients at 0 and (j)p^ we get 


f{Pl,P2)f{p2,P3)p2 - f{P2,Pl)f{pi,P3)pi + f{Pl,P3)f{P2,P3)p3 “ f {P2, P 3 ) f {f {Pl, P3)p3 + 

+9{Pl){f{P2,P3)) - 9{P2){f{PuP3)) + 

/(Pl.P3)(jl(P2)d^-^(/(p,.p.)|^ + 9(P2))(9.(p.)))=0. 

f{Pl,P2)9'{P2) - f{P2,Pl)9'{Pl) + [9{Pl),9{P2)] + 

(/(P2,Pl)^ + ^(P2)) (^l(Pl)))^(Pl) = 0. 

Expanding the first of these equations near the diagonal p 2 = P 3 and noting that 

f{Pl,P2)f{p2,P3)p2 + f{Pl,P3)f{P2,P3)p3 - f{P2,P3)f{PuP3)p3 = + 0 ( 1 ) 

P2 — P3 

we obtain 

9i{Pi)^^^^ -^(/(PnP2)^ + ^(pi))(^i(p2)) = 2 /(pi,P2)p2- 

diUid2Ui 9i[P2) V dp2 / 

Substituting this into our equations we arrive at relations (I3.19p . (Id.lSh for a local GT structure. 

One can check that all these steps are invertible and any local GT structure with relations 
fl3.18p . 03.191) gives a Gibbons-Tsarev system 05.35P with 

□ 

6 Integrability of Whitham type hierarchies 

In this Section we explain a relation between integrable Whitham type hierarchies and enhanced 
GT structures. 

Proposition 6.1. A Whitham type hierarchy with potentials hi{p,vi, ...,Vm), * = 1, ■■■,N 
is integrable by hydrodynamic reductions if and only if there exists a Gibbons-Tsarev system 
fl5.35p such that 

h'j{pi)di{hi{p2)) = h[{pi)di{hj{p2)), = ( 6 . 36 ) 
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by virtue of fl5.35p . 

Proof. The equation 


f{Pl,P2) = 


can be written as 

Er=i [KiPi)hj{p2)vi, - h'j{pi)hi{p2)y, 
h'Api)h[{p2) - h'Ap2)h[{pi) 


(6.37) 


and, therefore, coincides with the formula (77) from [19]. It is proven in [19] that the equation 
fl6.37p is equivalent to the integrability of a given Whitham type hierarchy. □ 

Proposition 6.2. There exists a one-to-one correspondence between integrable Whitham 
type hierarchies and enhanced local GT structures. Under this correspondence the space of 
potentials of a Whitham type hierarchy coincides with the space of solutions of the linear 
system fl3.26p . 

Proof. Write fl6.36p as ^ • ^7 executing di in numerators we get 

f{pi,P2)K{p2) + g{pi){hi{p2)) _ f{.Pl,P2)h'j{p2) + g{,Pi){,hj{,P2)) 


h'M 

Let A(pi,P 2 ) = 


h'Api 


, this function does not depend on i. Therefore, we get 


g{Pi){.hi{p2)) = A(pi,p2)h'(pi) - f{pi,P2)K{p2) 

which coincides with fl3.26p . Applying the relation fl3.18p to hi{p^) we can write 

g{pi)g{p2){hi{p^)) - g{p2)g{pi){hi{p:i)) = 

f{P2,Pi)h'{pi){hi{p:i))-f{pi,P2)h'{p2){hi{p^))+2f{p2,pi)pJ){pi){hi{p^))-2f{pi,p2)p^h{p2){hi{p:i)). 

Computing the l.h.s. and the r.h.s. of this relation by virtue of (I3.26p we obtain (I3.25p . □ 


7 The universal Whitham hierarchy 

In this Section we use notations introduced in Section 2, including G(p) and F{pi,p 2 ). 

According to [9] the universal Whitham hierarchy is given by potentials obtained by inte¬ 
gration of meromorohic differentials on a Riemann surface. We are going to construct such an 
hierarchy explicitly^ and prove that it is integrable by hydrodynamic reductions. 

Proposition 7.1. Fix constants Si,..., Sm such that Si-|-...-t-Sm = 1 (the simplest possibility 
is m = 1 and si = 1). The following formulas dehne an enhanced GT structure: 

” (9 d 

dip) = ^F{p,Uj)— + ^F{p,Wj )-— + G(p), f{pi,P2) = F{pi,p2), (7.38) 

j=i k=i 

^We need to choose constants of integrations carefully in order to obtain an integrable hierarchy. 
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Kpi,P2) 


E{pi,P2)p, 

E{pi,P2) 


E 

k=l 


Sk 


E{pi,Wk)pi 

E{pi,Wk) 


Moreover, the following functions belong to the space of potentials of this enhanced GT struc¬ 
ture: 

m 

hj{p)-hi{p), j = 2, Qaip) -^SkQaiwk), a = l,...,g 

k=l 


where 

m 

hj{p) = ln{E{p, Uj)) - ^ Sfc \-a{E{uj, Wk)). (7.39) 

k=l 


Proof. We need to prove identities fl3.25l) and (13.261) for given A(pi,p 2 ) and potentials. 
This can be done by straightforward computation using identities from Proposition 2.2. The 
simplest way is to start from identity fl3.26p for hj{p) — hi{p) and check it using identity (12.131) . 
It is clear from the proof of the Proposition 6.2 that (I3.25p is a consequence of (13.261) . It follows 
from Proposition 3.3 that ^ X{t,p)dt are also potentials of our hierarchy. Computing these 
integrals by virtue of (12.41) we conclude that the functions qa{p) — ^kQaiwk), a = 1, ...,g 
belong to the space of potentials. □ 

Proposition 7.2. The universal Whitham hierarchy is integrable by hydrodynamic reduc¬ 
tions. 

Proof. It is clear that the vector space spanned by derivatives with pespect to p of potentials 
described in the previous Proposition coincides with the space of meromorphic differentials on 
8 holomorphic outside Ui, ...,Un and with poles of order less or equal to one in these points. 
Therefore, we obtain a part of the universal Whitham hierarchy. In order to obtain the full 
hierarchy we apply the procedure of colliding point, see Proposition 3.2 and Remark 3.4. This 
proves the Proposition in the case g > 1. 

In the case = 0 we dehne an enhanced GT structure by (I3.29P and set X{pi,P 2 ) = 

The space of potentials contains the functions hj{p) — hi{p), j = 2, ...,n + 2 where hj{p) = 
ln(p — Uj), j = hn+i{p) = ln(p) and hn+ 2 {p) = ln(p — 1). This gives a part of the 

universal Whitham hierarchy corresponding to meromorphic differentials on CP^ with poles of 
order less or equal to one in ui, ..., 0,1. To obtain the full hierarchy we collide these points 

by a procedure similar to the one in the proof of Proposition 3.2, see also Remarks 3.3 and 3.4. 

In the case = 1 we dehne an enhanced GT structure by (13.301) and set 


A(pi,P 2) = p(pi -P2,r) -p(pi) - 2Tii. 

The space of potentials contains p — r and the functions hj{p) — hi{p), j = 2,...,n where 
hj{p) = ln(6*(p — Uj,T)) — \\i{6{uj,T)). This gives a part of the universal Whitham hierarchy 
corresponding to meromorphic differentials on 8 with poles of order less or equal to one in 
Ml, ...,Un- To obtain the full hierarchy we collide some of these points by a procedure similar 
to one in the proof of Proposition 3.2, see also Remark 3.4. □ 
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